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Abstract 

Angular distributions of fluorine and oxygen produced from 170 MeV 17 F 
incident on 208 Pb were measured. The elastic scattering data are in good 
agreement with optical model calculations using a double-folding potential 
and parameters similar to those obtained from 16 O+ 208 Pb. A large yield of 
oxygen was observed near #i a b = 36°. It is reproduced fairly well by a calcu- 
lation of the ( 17 F, 16 0) breakup, which is dominated by one-proton stripping 
reactions. The discrepancy between our previous coincidence measurement 
and theoretical predictions was resolved by including core absorption in the 
present calculation. 

PACS number(s): 25.60.-t, 25.60.Bx, 25.60.Gc, 25.70.-z 
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Our knowledge of atomic nuclei comes from experiments with nuclei in the valley of 
stability. The advent of radioactive beams enables us to explore nuclei near and beyond the 
drip lines and provides tests of current nuclear structure models. Nuclear reactions are often 
used as tools for studying nuclear structure. Breakup is an important reaction channel in 
the scattering of weakly bound nuclei and can be a rich source of information on reaction 
mechanisms and the structure of such nuclei [l|-[4j . 

Fluorine-17 is a proton drip-line nucleus with its valence proton bound by 0.6 MeV. 
Because of this loose binding, the r.m.s. radius, 3.7 fm ||, is significantly larger than 
that of the 16 core. Furthermore, the first excited state of 17 F, E x = 0.495 MeV and 
J 71 " = | + , is reported to have an extended r.m.s. radius, 5.3 fm |5[], and is considered to 
be a nuclear halo state. Studies of the influence on fusion of breakup of weakly bound 
nuclei has generated considerable interest in recent years [H]. The measurement of fusion 
of 17 F on 208 Pb at energies near the Coulomb barrier found no fusion enhancement Jr|. 
In our previous measurement of the breakup of 17 F by detecting the breakup products, 
proton and 16 0, simultaneously, a very small cross section was observed and disagreed with 
theoretical predictions ||. That experiment used a large area detector to optimize the 
collection efficiency; the measured cross sections were averaged over the angles spanned by 
the detector. 

A recent study of the breakup reaction 8 B — > 7 Be+p on a 58 Ni target at sub-Coulomb 
energies [[| demonstrated that the resulting 7 Be angular distribution can be accounted for by 
two different theoretical approaches, coupled discretized-continuum channels (CDCC) Jirj 



and dynamical calculation |fLl|| . With this in mind, we have performed a new experiment 
on the breakup of 17 F by measuring the angular distribution of reaction products in singles. 
The data provide some details of the cross section over a range of angles and no Monte Carlo 
simulation is required for obtaining the efficiency of the detectors, unlike the coincidence 
measurement. However, the identification of reaction products with the breakup channel 
is not as straightforward as in the coincidence measurement, since similar products can be 
produced in several reaction channels. With the aid of theoretical calculations, the dominant 
reaction products were shown to originate from breakup. 

The experiment was carried out at the Holifield Radioactive Ion Beam Facility (HRIBF) 
using a 170 MeV 17 F beam incident on a 2 mg/cm 2 208 Pb target. A 44 MeV deuteron beam 
provided by the Oak Ridge Isochronous Cyclotron (ORIC) was used to bombard a fibrous 
hafnium oxide target to produce radioactive 17 F by the 16 0(d,n) reaction |12|]. The A=17 



ions were extracted from the target-ion source, mass analyzed and subsequently accelerated 
by the 25 MV tandem postaccelerator. The 17 isobar contaminant was removed by inserting 
a carbon stripper foil at the exit of the tandem accelerator, before the 90° analyzing magnet, 
and selecting the fully stripped 17 F 9+ ion beam. The beam was monitored by a Si surface 
barrier detector positioned at 10° to the beam direction in the target chamber as well as by 
the focal plane detector in the Enge spectrograph at 3° on the other side of the beam. The 
average beam intensity was 5 x 10 5 17 F/s. The reaction products were detected at 9^ = 45° 
by a AE-E telescope consisting of a large area (900 mm 2 ) 100 fim Si surface barrier detector 
(SBD) and a 5 cmx5 cm, 300 fim double-sided Si strip detector (DSSD) behind the SBD. 
The DSSD, which has 16 vertical and 16 horizontal strips, was placed 8.3 cm from the target, 
resulting in an angular resolution of ~ 2°. At forward angles, the elastically scattered 17 F 
was measured by the DSSD only. 
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Figure [T] shows the E vs. AE histogram of the SBD-DSSD telescope for events in all the 
strips. Reaction products of Z=9(F), 8(0), 7(N), 6(C), 5(B) and 1(H) were observed. The 
band of constant energy loss (AE) at channel number ~ 135, corresponding to the energy 
loss of the elastically scattered 17 F, results from the positron decay of 17 F stopped in the 
DSSD. These events were recorded when a positron emission took place at the same time as 
an elastic 17 F struck another strip in the DSSD. 

The elastic scattering data were extracted from the Z=9 products identified in this E 
vs. AE histogram. The energy resolution of the detectors was not good enough to allow 
a clear separation between the elastic and inelastic scattering events. Moreover, the mass 
of the reaction products cannot be identified in this experiment. Since there are some 
neutron pickup channels with positive Q values, the elastic scattering data may include 
contributions from inelastic scattering and neutron transfer reactions. The absolute cross 
section was obtained by normalizing the yields to the Si detector at 10° where the elastic 
scattering was taken as Rutherford scattering. The elastic scattering angular distribution 
is shown in Fig. 0. As can be seen, the assumption of pure Rutherford scattering at 10° is 
valid. 

Optical model fits to the elastic scattering data were performed using the code Ptolemy 
T3fl . The Woods-Saxon potential parameters were taken from the 192 MeV 16 + 208 Pb 



elastic scattering [14] and only the depth of the imaginary potential was allowed to vary. The 
best-fit result is shown by the dotted curve in Fig. § and the parameters are listed in Table | 
as SET I. Parameters in SET II of Table | were obtained using the procedures described in 
Ref. |14| by fitting the data with the depth of the real potential, V, fixed at 40 MeV and 



varying the other parameters. Throughout the fitting processes, the shape of the real and 
imaginary potential was assumed to be the same. The imaginary potential depth, W, was 
changed from 35 to 95 MeV in 5 MeV steps and the radius and diffuseness parameters, r 
and a, were varied to find the minimum \ 2 . The best fit is shown by the dash-dotted curve. 
The radius and diffuseness parameters were very similar to those obtained in the 140 MeV 
i6q _j_ 208pk e i as ti c scattering [Tj|. The dashed curve, corresponding to parameter SET III, 
is the result of allowing all the parameters, V, W, r, and a, to vary subject to the constraint 
that r and a were identical for the real and imaginary potentials. The dashed curve and 
dash-dotted curve are almost indistinguishable. All three sets of parameters describe the 
data very well at angles smaller than 41°. However, differences between the dotted curve 
and dashed curve can be seen at larger angles. The inclusion of inelastic scattering and 
transfer in the data may account for the fact that the experimental points are above the 
dotted prediction. 

The elastic scattering data were also compared to calculations using potentials obtained 
from a double-folding model. The effective nucleon-nucleon interaction was taken from 
Satchler's systematics [|16|] which uses the three-parameter Yukawa form 



U NN = -(V + iW)^ 

s/t 

where s is the distance between the two nucleons and t = 0.7 fm. The depth of the potentials 
was V = W = 60 — 0.3E/A = 57 MeV where E/A is the energy per nucleon of the projectile. 
The nuclear density distribution of 208 Pb was constructed using the two-parameter Fermi 
form with parameters given in Table 2 of Ref. [TE\. The ground state density distribution 
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of 17 F was constructed from shell model single particle wave functions as described by 
Satchler [|17|] (method B). Calculations using this double-folding potential with no adjustable 
parameters are in very good agreement with the data, as shown by the solid curve in Fig. |3|. 
Beyond about 45°, the calculation underpredicts the data somewhat, which may be due to 
the inclusion of inelastic scattering and neutron transfer reactions in the data. The inelastic 
excitation of 17 F to its first excited state, E x = 0.495 MeV and J n = | + , and 208 Pb to its 
first 2+ E x = 4.09 MeV, and 3~, E x = 2.61 MeV, states were calculated in the DWBA by 
the code Ptolemy. The sum of the inelastic scattering predicted by DWBA and the elastic 
scattering predicted by the folding potential is shown by the dotted curve in Fig. [5] which is in 
excellent agreement with the data. Among the three inelastic channels, the excitation cross 
section for populating the first excited state of 17 F is calculated to be the largest, followed 
by exciting 208 Pb to its first 3~ state. Since the single particle wave function used seems 
to give a good account of the ground state density distribution of 17 F, the halo property of 
the first excited state in 17 F could be studied if the inelastic scattering could be resolved in 
the data. Coupled-channels calculations were performed with these three inelastic channels 
using Ptolemy. As can be seen by the dashed curve in Fig. ^ the elastic scattering is 
influenced only slightly by the channel couplings. The one-neutron transfer reaction was 
calculated by DWBA using Ptolemy. The transfer yield was predicted to be at least two 
orders of magnitude less than the elastic scattering and can be ignored, consistent with the 
agreement between the data and the dotted curve in Fig. ||. 

The angular distribution of Z=8 reaction products used the same normalization as the 
elastic scattering. Since this experiment provided no mass identification of the reaction 
products, the observed oxygen can originate from nucleon transfer, charge exchange, and 
breakup. The one-proton transfer has the largest Q value of all transfer channels lead- 
ing to an oxygen isotope, 16 0, as the final product. Finite range DWBA calculations us- 
ing the code Ptolemy were performed to predict the cross section angular distribution of 
208 Pb( 17 F, 16 O) 209 Bi. Transfer to the lowest six single particle states in 209 Bi was calculated. 
The spectroscopic factors were assumed to be 1 which will give upper limits for the cross 
sections. The sum of the six cross sections is shown by the dotted curve in Fig. |j. It can be 
seen that the one-proton transfer cannot account for the yields of the oxygen angular distri- 
bution. Similar results were obtained for calculations using the code FRESCO ||18|| . Although 
reactions leading to 18 and 207 Bi in the exit channel have positive Q values, they cannot 
occur by simple single-step transfer processes. Therefore, the cross sections are expected to 
be smaller than that of one-proton transfer |T9 . 



Oxygen-17 can be produced by the charge exchange reaction, 208 Pb( 17 F, 17 O) 208 Bi, with 
a Q value of -0.11 MeV. Calculations were performed to estimate the cross section for this 
reaction channel. Charge exchange reactions can take place by direct charge exchange, a 
combination of inelastic excitation and direct charge exchange, or by successive one-nucleon 
exchange. Since the direct charge exchange is orders of magnitude smaller than one-nucleon 



transfer ||20|| , calculations were carried out using FRESCO to model the successive one-nucleon 
exchange processes [2"T]. Two possible modes, proton stripping followed by neutron pickup 
and neutron pickup followed by proton stripping, were considered. Only transfer between 
ground states was calculated. The charge exchange cross section is calculated to be 2 orders 
of magnitude smaller than that of one-proton transfer calculated by DWBA. Measurements 
of 28 Si( 18 0, 18 F) 28 Al show that the charge exchange cross section is at least an order of 
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magnitude smaller than that of one-nucleon transfer |22| , consistent with the results of our 
calculations. Since DWBA calculations predict that the one-proton transfer cross sections 
are much smaller than the measured cross sections, it is safe to assume that charge exchange 
can be ignored. 

The measurement is inclusive, so products from two breakup processes, diffraction and 
stripping, can contribute to the data @. In diffraction dissociation, the projectile breaks 
up, leaving the valence nucleon in the continuum and the core intact. The final state of the 
reaction consists of the valence nucleon, the core of the projectile and the target nucleus in 
its ground state. This is the process that our previous coincidence experiment measured || . 
The stripping breakup depends on the core-target and nucleon-target absorption potentials. 
In a loosely bound nuclear system, if the separation energy (Sn) is much less than the 
kinetic energy per nucleon, Sn << E/A, the valence nucleon and the core can be treated 
approximately as independent particles p3fl . The effects of nucleon-target and core-target 



interactions become noticeable if the interactions are strong. 

Semiclassical calculations of breakup similar to those in Ref. |2| were performed to inves- 
tigate these processes. The wave function of the projectile after interacting with the target 
is described by the ground state wave function of the projectile and the profile functions of 
the proton and 16 core. The interaction between the target nucleus and the constituents of 
the projectile was obtained from Ref. [24]. The 16 0-target optical potential was calculated 



in a single-folding model based on the CH89 optical model potential ||24|| . The breakup 
probability as a function of the impact parameter was calculated and converted to an angu- 
lar distribution assuming pure Coulomb scattering. The calculated angular distribution was 
converted from the center of mass of the reaction to the laboratory frame using the Jacobian 
for elastic scattering. Because the reaction energy is high, the uncertainty introduced by 
this approximation is smaller than the angular resolution of the DSSD. 

Figure f| presents the comparison of the experimental oxygen cross section data with 
model calculations. The long-dashed curve is the prediction for stripping breakup while the 
short-dashed curve is for diffraction breakup. Apparently, stripping breakup is the domi- 
nant reaction in the measured angular range. The solid curve shows the sum of diffraction 
and stripping. It overpredicts the measured cross sections slightly; however, the data do 
not extend far enough forward for a good comparison. It should be pointed out that the 
calculations employed here did not consider the recoil of the 16 0. In addition, excitation of 
the 16 and 208 Pb in the final state was ignored. Measurements of the breakup of 17 C by 
197 Au show that core excitation is important for nuclear breakup pSfl , but very small for 
Coulomb breakup pTJ. In our case, the measurement was performed near the grazing angle 
where the nuclear breakup should be significant. Therefore, the excitation of the core and 
target nuclei should be considered. Furthermore, the stripping and diffraction breakup were 
calculated separately. A more complete theoretical treatment of this subject is underway 
and will be published elsewhere ||27|| . 

Diffraction breakup of 17 F was measured previously in our coincidence experiment ||. A 
first-order perturbation calculation overpredicted the measured cross section by a factor of 
4. The present work shows the importance of stripping breakup resulting from the nucleon- 
target interaction. The discrepancy in the coincidence experiment can be accounted for by 
including the core absorption by the target nucleus which was not considered in the previous 
calculations. As shown in Fig. [5], the data agree with the calculations taking into account 
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the core-target interaction (short-dashed curve). The calculation is also consistent with 
predictions made by the CDCC calculation [|l(]] shown by the dash-dotted curve. The CDCC 
calculation takes into account of higher order effects, continuum-continuum couplings, and 
full three-body kinematics to produce results in the laboratory frame |28 . 



It is shown in Ref. [Q] that the breakup of n Be (Sn = 0.5 MeV) by 208 Pb is dominated by 
stripping at 800 MeV/nucleon. At 40 MeV/nucleon, the stripping and diffraction breakup 
cross sections become similar 0. In the breakup of 8 B at energies below the Coulomb 
barrier, the stripping and diffraction breakup cross sections are also similar ||. In this 
work, it is found that the stripping breakup cross section of 10 MeV/nucleon 17 F is almost 
twice as large as that of the diffraction breakup. 

In summary, the elastic scattering and breakup of 17 F was measured on a 208 Pb target. 
The elastic scattering data were fitted with an optical model using a Woods-Saxon potential. 
The potential parameters are very similar to those obtained from 16 O+ 208 Pb elastic scatter- 
ing. The potential generated from a double-folding model using Satchler's systematics also 
reproduced the elastic scattering data. The dominant contribution to the measured angular 
distribution of oxygen nuclei in the exit channel was found to be the one-proton stripping 
reaction resulting in 16 0. This demonstrates the importance of the proton-target interaction 
in the breakup of 17 F on 208 Pb. Based on this finding, calculations of diffraction breakup 
considering core absorption were carried out. The new calculations agree with our previous 
coincidence measurement. 
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TABLES 





V 


W 


r = ri 


3i — Sii 




(MeV) 


(MeV) 


(fm) 


(fm) 


SET I 


40 


49 


1.23 


0.63 


SET II 


40 


75 


1.28 


0.50 


SET III 


50 


93 


1.27 


0.51 



TABLE I. Optical model potential parameters obtained from fitting the 170 MeV 17 F+ 208 Pb 
elastic scattering data. The depth of the real and imaginary potentials is shown by V and W, 
respectively. The real and imaginary radius parameters are shown by r and r;, respectively. The 
real and imaginary diffuseness parameters are shown by a and a;, respectively. 
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FIGURES 



FIG. 1. Histogram of E vs. AE obtained by the SBD-DSSD telescope. 

FIG. 2. Angular distribution of the ratio of measured elastic scattering to calculated Rutherford 
scattering. The optical potential fit to the data is shown by the dotted curve for parameters SET I, 
dash-dotted curve for SET II, and dashed curve for SET III. The dash-dotted curve and dashed 
curve are almost indistinguishable. 

FIG. 3. Angular distribution of the ratio of measured elastic scattering to calculated Rutherford 
scattering. The solid curve is the result of calculations using the optical potential obtained from 
a double-folding model. The sum of elastic and inelastic scattering calculations is shown by the 
dotted curve. The result of coupled-channels calculations is shown by the dashed curve. 

FIG. 4. Angular distribution of oxygen. The long-dashed and short-dashed curves are for strip- 
ping and diffraction breakup, respectively. The dotted curve is the result of a one-proton transfer 
DWBA calculation. The solid curve represents the sum of stripping and diffraction breakup. 

FIG. 5. Data and predictions of diffraction breakup. The long-dashed curve is for the calcu- 
lation without considering core absorption and the short-dashed curve includes core absorption. 
The CDCC calculation is shown by the dash-dotted curve. The open point is from our previous 
measurement ||. 
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